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We investigate (non-relativistic) atomic systems interacting with quantum electromagnetic field (QEF). The 
resulting model describes spontaneous emission of light from a two-level atom surrounded by various initial 
states of the QEF. We assume that the quantum field interacts with the atom via the standard, minimal- 
coupling Hamiltonian, with the A 2 term neglected. We also assume that there will appear at most single 
excitations (photons). By conducting the analysis on a general level we allow for an arbitrary initial state 
of the QEF (which can be for instance: the vacuum, the ground state in a cavity, or the squeezed state). We 
derive a Volterra-type equation which governs the time evolution of the amplitude of the excited state. The 
two-point function of the initial state of the QEF, integrated with a combination of atomic wavefunctions, 
forms the kernel of this equation. 
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1 Introduction 



The phenomenon of emission of photons during the transition between atomic energy-levels is one of the 
most important, and well understood aspect of contemporary theoretical physics. The explanations of this 
process which have already been presented in the early days of quantum theory [1], vary in their level 
of sophistication. The desirable formulation, in the spirit of quantum theory, describes the process as 
a unitary evolution of the full system which is composed of the atom and the quantum electromagnetic 
field. It is disappointing that modern literature in most cases reports models which either do not fit to this 
basic requirement of quantum theory, or at some point involve unjustified operations with (unnecessarily) 
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singular expressions 1 . One has an unpleasant feeling that the outcome (the exponential decay of the excited- 
state's amplitude) is simply assumed at some point of the investigations. While we do not claim that no 
clean model exists, we believe a simple field-theoretical model which is additionally more general with 
respect to the quantum electromagnetic field (QEF), is worth presenting. 

It is generally recognized that the atomic emission of radiation is influenced by both: the atomic struc- 
ture (in particular: the energies and wavefunctions of the levels under consideration) and the initial state of 
the QEF (which can be a strong, coherent state or a vacuum, for instance). 

With respect to the latter we note that there exists a great variety of states of QEF. On the level of local 
quantum field theory [3], it is known that even disjoint representations of the algebra of fields (the Canon- 
ical Commutation Relations (CCR)-algebra) can be constructed 2 . Each such representation carries its own 
notion of excitations ("photons"). A general method of Gelfand, Naimark and Segal (GNS) allows to con- 
struct a representation once a single "reference state" £1 is specified. If we deal with fields propagating on 
the Minkowski spacetime, the Fock representation, based on the usual vacuum, is physically distinguished. 
However, if we were interested in some generic spacetimes (or time-dependent environments) there would 
be no privileged "reference state". On the mathematical level, one is content with a picture of many per- 
missible representations, each based on a different "reference state" f2, although such a state of affairs is 
certainly conceptually worrying (there is no unique answer to the question what a photon is). 

On the other hand, experimental advances in the area of quantum optics have produced a variety of 
states of light which exhibit non-classical features. The so-called squeezed states or the ground states 
of the QEF in small cavities provide environments which are much different form the usual Minkowski- 
spacetime vacuum. In particular, the experiments with atoms surrounded by squeezed states, [4], revealed 
a modification of their spectroscopic characteristics. 

One of the simplest ways to test the properties of various states of QEF is to let the quantum field in- 
teract with a quantum system of relatively simple structure. As an example, we can consider a "two-level 
atom", assume a simple minimal-coupling interaction with the QEF, and investigate the time-evolution of 
various initial states of the full system. Such an approach is also experimentally viable, as this evolution 
is nowadays fully accessible to measurements (see e.g. [5]). In order to illustrate the potential of such an 
approach we remark on a possibility to measure the Unruh effect with the help of an electron, namely, a 
simple perturbative calculation shows that the state of a two-level quantum system coupled to the QEF ini- 
tially in the KMS (thermal) state thermalizes to the temperature of the KMS state. As originally suggested 
by Bell and Leinaas [6], a spin of a linearly-accelerated electron might serve as a detector of the Unruh 

1 See, for instance, [2] chapter 6.3. 

2 It is possible, for instance, to prescribe two states of the QEF which cannot be expressed as density-operator states in a common 
Fock space. 
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temperature 3 . It is one of the aims of this paper to pursue further this type of models in which test atomic 
systems are employed as detectors of properties of various states of the QEF. 

The paper contains a field-theoretical model of an interaction of a two-level atom with the QEF. We 
assume the atom to be initially in an excited state, and let it interact with the QEF; the dynamics of such 
a system is usually referred to as the spontaneous emission of radiation. It has been investigated in the 
past by many authors with various emphasis on mathematical rigor, on the one hand, and on physical 
concreteness on the other. We shall briefly summarize these attempts here, in order to put our model into a 
context. 

The simplest way to estimate the dynamics of the full system is to use the time-dependent perturbation 
theory w.r.t. the interaction of the QEF and the atom (Eq. (3)). By doing so, one can approximate the initial 
dynamics of the emission of a photon, however, the long-term behavior of this process cannot be estimated 
in this way. Another method, due to Weisskopf and Wigner [l] 4 , establishes an ordinary integro-differential 
equation for the time-dependence of the amplitude of the excited state. The model we shall present here 
generalizes the Weisskopf- Wigner model with respect to the allowed initial states of the quantum radiation 
field. 

With an alternative approach, Bach and collaborators [7] have investigated systems of non-relativistic 
particles coupled to the QEF from the functional-analytic point of view. Without any assumption on the 
number of atomic energy levels, they worked with the full minimal coupling Hamiltonian (together with 
the quadratic term, and an ultraviolet cutoff), and have been able to prove the existence of a ground state 
of the coupled system. Moreover, an estimate has been derived for the life-times of the states of the atom 
which were stationary when the interaction was absent. 

A common feature of all these developments is, that the choice of the initial state of the QEF is made at 
a very early stage of investigations. Consequently, given the final equation for the excited state's amplitude, 
it is difficult to discern which features of this equation are related to the initial state of the QEF, and which 
are influenced by the atomic wavefunctions 5 . 

In this paper we attempt to generalize the Weisskopf- Wigner model; in particular we proceed in a 
general way, so that the initial state of the QEF can be left unspecified to the very end. By doing so, we 
prove that the spontaneous emission process is governed by a Volterra integro-differential equation (Eq. 

The usual vacuum state, when observed by an accelerated observer, appears to be a thermal state, with an acceleration- 
dependent temperature. 

4 See, for instance, the textbook of Scully and Zubairy [2]. 

5 In an interesting paper Bondarev and Lambin [8], investigate the spontaneous emission of light from atoms placed inside of 
small cavities (carbon nanotubes). Their explicit numerical results complement our analysis and exemplify the possible complexity 
of the emission process. 
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(15)), the kernel of which is composed of the two-point function, (fi, Ai(x)Aj (y)ti), of the initial state of 
the radiation field, smeared with test functions which are derived from the wavefunctions of the ground- 
and excited-state of the atom. As the two-point functions are the foundation of the GNS construction, 
and on the other hand are known for physically interesting states, this opens up a possibility to investigate 
various initial states of radiation on an equal basis. 



2 Formulation of the model 

2.1 General setting 

A natural Hilbert space for a system of an atom (non-relativistic) coupled to the quantum radiation field is 

W = L 2 (R 3 ) <g) JT, (1) 

where T denotes the physical (transversal) Fock space of the radiation field. This space will be later 
restricted to a subspace, consisting of a two-dimensional Hilbert space of atomic states, tensor-multiplied 
with the vacuum- and one-excitation-subspace of T . The single excitations (over the reference state £1) of 
the quantum radiation field, denoted by the symbol | /), will be described by 



|/) = A(f )n = J d A x / j (x)A*(0,x)fi, (2) 

where ^4,(0, x) is the (unsmeared) vector-potential operator at t = 0. The complex-valued test function 
/i(x) will carry the information about the "wave packet" of the excitation 6 . 
The interaction between the atom and the QEF will be generated by 

V(t,x) = -V^(t,x)- Pj , (3) 

where 7 pi = —idi, a is the fine-structure constant, and A z (t, x) denotes the vector-potential operator, 
in the radiation gauge, with the free time-evolution already implemented 8 . This is a standard interaction 
Hamiltonian with the A 2 term neglected. 



6 If ri is the vacuum on the Minkowski spacetime, the excitations are photons. In this case A(f)£l = J f (p) a a(p) ^> 
where a* (p) is the creation operator for o-th polarization and f a (p) = J d 3 x p (x)ef (p)e _lpx . The symbols ef (p) denote the 
transversal polarization vectors. 

7 In the paper we use dimensionless units (see appendix A). By the letter p we denote the frequency: p = |p|. 

8 Which means that the free time-evolution of A' (t, x) is governed by the Maxwell equations. 
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2.2 Dynamics of the restricted system 

In this section we will consider the evolution of the above system in the situation where the radiation field 
is initially in the state ft, about which we only assume that the expectation value of an odd number of field 
operators in this state vanishes 9 , and the atom is initially excited. 

The HUbert space of the atomic motion will be restricted to a two-dimensional subspace (two energy 
levels), and the Fock space will be restricted to the zero- and one-excitation-subspace 10 . 

We start by restricting the interaction (3) to the subspace of only two energy-levels. To do so, we 
substitute the interaction operator V by V 2 : 

v 2 = Wi'Vilij)* iv^XM 

y=0,l 

{tpi,V i)j) L 2 = J d 3 x Vi(x)^(x), 
where we have employed the Dirac's notation. We obtain 11 : 

v 2 = i^[\i> 1 )(M®Mx) + hto><iM®^t fa)], (4) 

where 

Xi(x) = Vi(x) <9jVo(x), 
ipi(x) = V'o(x) ^Vi(x) 

Mx) = J d 3 xMt,x) X l (x). 

Because of A t (x) = —A t (ip), the interaction can be written as 

V 2 = Va<T 2 ®A t (x), (5) 

which is a selfadjoint operator on C 2 ® T . At this point we remark that the process of "spontaneous 
emission of radiation" will be generated by a selfadjoint Hamiltonian, and thus will be governed by a 
unitary evolution. In this way the basic requirement of the quantum theory will be fulfilled. 

We shall now additionally restrict the investigations to the subspace of at most one excitation of the 
QEF 12 . Later, we will employ the interaction picture and derive a system of evolution equations. At the 



9 The vacuum, the ground state, the squeezed state and a general class of quasi-free states fulfill this requirement. 

10 Every vector v of such a space can be written as: v = cfl + A(f)U with a complex number c and complex, vector-valued 
functions /*(x). 

11 The terms and |i/'o}(V'o| drop out, if the wavefunctions are real (which we hereby assume), because 
/ 4>o(di'4>o)A'' = 0, as a consequence of the chosen gauge condition diA 1 = 0. 

12 This restriction allows us to derive a single equation for the evolution of the excited-state's amplitude (Eq. (15)). We expect 
such a restriction to be justified in some cases, e.g. where one-photon transitions are known (from experiments) to be dominant. 
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time t, the state of the full system, in the restricted space, can always be parameterized by two complex, 
vector- valued functions f t (x) , g t (x) (excitation wavepackets), and two complex numbers c(t), d(t): 

S = c(t) tpx <g> Q + d(t) O ^ + tpo ® A (f t )Q + Vi ® A (g t )Q. (6) 

We note that we speak of excitations because the state SI is still unspecified. We use the following initial 
conditions 

c(0) = 1, 
/o(x),so(x),d(0)=0. 

For such conditions the functions g t (x) and d(t) never acquire non-zero values, and will be omitted from 
now on. In the interaction picture, we reparameterize the state of the full system as follows: 

Si(t)=c(i)^i®fi + Vo®A(/t)fi. (7) 

The time evolution of c(t) and f t needs to be determined; it will be generated by the operator V 2 7 , which is 
the operator V2 in the interaction picture: 

Vl = 1^ (e^lViXVol - e-^IVoXV-il) ® ^(X), (8) 

where a; = £q — .En is the energy difference of both states of the atom. We shall introduce a decomposition 

of ft in terms of an orthonormal basis {/,„} of the one-excitation space 

00 

A(f t )n = c m (t)A(f m )n, (9) 

m—l 

where f m are complex functions, orthonormal w.r.t. the scalar product 13 

(f,g) = f A Wl3 (0,x,0,y)f (xV(y). (10) 



The evolution equation, 

i^-(t) = V 2 I (t)S I (t), (11) 

after taking scalar products with ^i(g)0 and ip ®A(f n )Sl, leads to an infinite system of ordinary differential 
equations: 

00 

ic(t) = e iojt ■ c m (t) (n,A t ( x )A(f m )n), (12a) 

m—l 

ic n (t) = -1^ e"*" ■ c(t) (SI, A(f n )A t ( X )n). (12b) 



13 Here, the two-point function of the initial state of the QEF, ujij (t, x, s, y) = (f!, Ai (t , x) Aj (s, y)Q) , enters the investiga- 
tions. The one-excitation space is separable for reasonable initial states fi. 
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Integrating the second equation(s), (12b), from to t, with the initial condition c„(0) = (for all n), yields 

c n (t) = -yfc f ds e-^ s c(s) (O, A(f n )A s ( X )n). (13) 
Jo 

After substituting this into the equation for c(t), (12a), we obtain 

,.t oo 

c(t) = -a ds c(s) e^-^ V (fi, A(f m )A s ( X )n) (O, ^ t ( X )^(/ ro )0), (14) 

m=l 

and finally, as a consequence of the completeness of the {/ m } in the one-excitation subspace of T, we get 

c(t) = -a[ dsc{s)e l ^ t - s \n,A t (x)A s ( X )n) = -a [ ds c{s)e iuj( - t - s '> S(t, s). (15) 
Jo Jo 

with 

S(t,s) = J d 3 xd 3 ycu lJ (t, Xl s 1 y) X l ^)x l (y)- (16) 

The equation (15), valid for arbitrary initial states of the QEF, is a Volterra integro-differential equation 
of the second type. In the important case of a stationary initial state ft, the smeared two-point function 
s 7 x) depends only on the difference 14 , t — s, and therefore in this case we obtain a convolution 
Volterra integro-differential equation. Such an equation can be brought into the form of an integral equa- 
tion. To show this, we integrate it w.r.t. the variable t from zero to T and (with an appropriate change of 
variables) obtain 

c(T) = 1 - / Z(T - s)c(s) ds, (17) 

where 

Z(t) = a [ dt S{t,ti)e iut . (18) 



An alternative way to approach the equation (15), in the case of stationary initial states, is to perform a 
Laplace transform of it. For the transformed amplitude, c(s), we obtain a relation 

c{s) = — , (19) 

s + ab{s — iuu) 

where S(s) = f °° exp(— st) S(t, 0) dt. The inverse Laplace transform gives the original amplitude which 
we intend to compute: 

1 f e st 

c(t) = — / ds, (20) 

2?r« Jc s + aS(s - iuj) 



14 As a consequence, S(t, s) = S(t — s, 0). 
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where the contour of integration, C, needs to be chosen parallel to the imaginary axis, to the right of all 
poles. In our case, the dynamics of the total system is unitary, and consequently there will be no poles 
of the integrand for Re(s) > 0. The standard procedure is to attempt to close the contour C from the 
left (Re(s) < 0). The dominant long-term asymptotic behavior of c(t) comes from the singularities of 
c(s) with the greatest Re(s). In general, the integrand, analytically continued to Re(s) < 0, may exhibit 
branch points (which often lead to power-law decays) and simple poles (which lead to exponential decays). 
It is not possible to be any more specific here, without specifying ft, as the analytic properties S(s — iuo) 
obviously depend crucially on it. 

In any case, even for non-stationary initial states, the kernel of the equation (15) is a smooth function 
of its arguments (because the test functions \ are smooth), and therefore our problem always reduces 
to solving a Volterra-type integro-differential equation of the second kind with a smooth kernel. The 
numerical experience with this equation, for various initial states, shows that the solutions exhibit a non- 
trivial initial dynamics which is followed by a relatively uncomplicated, monotonic-decay phase 15 . 

We wish to defer a systematic study of the solutions of (15) to a future publication [11]. 

2.3 Special case: Minkowski vacuum as the initial state ft 

The simplest case of a great physical importance arises, if one takes the usual vacuum state £1 as the 
initial state of QEF. Then, the electromagnetic field operator can be expressed in terms of the creation and 
annihilation operators, 

d 3 k 

{Zi) 

2tT 



1 f d A k 

Mt, X ) = J — e f (k) {a* (k) <? kx + a a (k) e~ ik *} , 



where ef (k) denote the two (a = 1,2) k-dependent, transversal polarization vectors, and the vacuum is 
defined via 

a a (k)n = 0. (22) 
The respective two-point function is well-known: 

w?-(t,x, s , y ) = (n,Mt,x)M*,y) °) = / $ e ^ s) (*« - ^r) e4p(x_y) - (23) 

Denoting by Xi(p) tne Fourier transform of Xi( x )> 

Xi(p)= /Vze-^^x^W, (24) 

15 If the atom is placed inside of a reflecting cavity, however, c(t) exhibits a completely different dynamics, due to the interaction 
with the emitted photon. In an extreme case this can even lead to oscillations of |c(t)| 2 which are called the vacuum Rabi oscillations 
[8, 9]. 
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we get 

S(t ~s) = ^ / ^ e-«*-) (V - ^(p) (25) 

The function — s) is evidently a smooth function of its arguments 16 , because x(p) is a smooth function 
of rapid decay 17 . 



2.4 Another special case: squeezed state as the initial state £1 

The squeezed states are examples of states of the QEF which can be produced experimentally [2], and 
exhibit exotic, inherently quantum properties 18 . These states appear as a coherent superposition of pairs 
of photons. The photons constituting a pair, in the case of a degenerate squeezed state, possess the same 
wavepacket, f a (p). (We assume this function to be smooth and rapidly decaying.) The squeezed states are 
important in the context of this paper, as their presence is known (experimentally) to influence the atomic 
dynamics [4]. 

A squeezed state, \S), is a vector state created in the Fock space constructed upon the (distinguished, 
Minkowski-spacetime) vacuum with the help of a unitary operator S(f, r): 

\S) = S(f,r)Q = e r [ a W') a («" a *W) a *W)]f7, (26) 

where r is a real number which we interpret as the amplitude of squeezed light 19 . The operator a(f) 
denotes here the smeared annihilation operator 20 : 

«(/) = f d 3 pjjp)a a (p). (27) 

The following (commutation) relations hold for an arbitrary test function g(p), and facilitate the calculation 
of various expectation values: 

S*a(g)S = a(g) - (g, f) a(f) + (g, f) [a(f) cosh r - a* (/) sinh r] (28) 
S*a*(g)S = a* (g) - (f,g) a(f) + (f,g) [a*(/)coshr - a(/)sinhr] , (29) 



16 In particular, the limit s — > t is finite. 

17 The smoothness of x(p) assures, that the point p = is an integrable singularity. The decay property for large |p| makes the 
integral convergent for large momenta, even when the time derivatives of S(t — s) are taken. 

Such as, for instance, a sub-vacuum level of electric field fluctuations. 

19 As the arguments of S(f, r) will not be varied in the following discussion, we will briefly write S instead of £(/, r). 
20 

The index a is summed over, a = 1,2, whenever it appears twice. 
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where (/, g) = J d 3 p f a (p)g a {p) stands for the scalar product of f a and g a . With the help of these 
relations we find 



(Sn,a* a (p)a (k) SO) = f a (p)U(k) sinh 2 r 
(SQ,a a (p)ap(k) Sfl) = -/ a (p)/a(p) sinhrcoshr 

The two-point function of the radiation field, in the squeezed state, can now be computed (here x = (x , x)) 

cjf j (x,y) = (Q,A i (x)A j (y) Q)=^ j (x,y) + J ^0 e «(p)ef (k)- 
• {-[f a {p)fp(K)e- ipx - iky +c.c] sinhrcoshr + \Mp)U(k)e ipx - iky + c.c] sinh 2 r} . 

Evidently, it differs from the vacuum-two-point function, oj^(x, y), by a function which is smooth in both 
arguments. 

If the wavepacket f a (p) is concentrated around a specific momentum q, then the two-point function of 
the squeezed state simplifies to 

ujfj(x,y) =uj° J (x,y)-{\d l d :j e- iq(x+y '> + c.c. sinhrcoshr + T^e 1 ^-^ + c.c. sinh 2 r}, 

where d is a polarization vector (orthogonal to q). From the above expression one realizes that the dif- 
ference between lu s and uj° becomes negative for some periods of time. This fact corresponds to the 
sub-vacuum fluctuations of the vector potential 21 Ai\ the expectation values of the squared electric-field 
operator, and the energy-density operator behave similarly. We note that the squeezed state is not a sta- 
tionary state, and therefore the Volterra equation (15) is will not be a convolution-type equation in this 
case. 

3 Example 

3.1 Hydrogen atom, 2P — > IS* transition 

The wavepackets of the initial and final state of the atom are 

a 3 / 2 a 5 / 2 

4V27T 



We will compute the smearing functions, 

Xi(p) = J d 3 xe-*P* ^(x) a^o(x). 
21 The square of A l (x) is defined via normal ordering w.r.t. the vacuum fi. 



(30) 
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The differentiation w.r.t. x l leads to a factor x l (multiplied with functions of r only) which can be expressed 
as a derivative id/dpi. Similarly, the term z — rcos(9), can be expressed as a derivative w.r.t. the p z . In 
this way we obtain 



Xi(p) 



a" 



d 2 



4V2tt dp z dp,, 



I 



d 3 xe~ ipx - 



-3/2ar 



and consequently 

Xi (p) = a 5 V2 



4 Pi p z 



(31) 



(32) 



_(p 2 +9/4a 2 ) 3 (p 2 + 9/4a 2 ) 2 _ 

The first term in the square bracket is longitudinal (proportional to pi) and will be canceled by the projection 
on the transversal directions in the two-point function. 



3.2 Spontaneous emission in the presence of vacuum 

We now compute the kernel of the equation (15) for the 2P — > 15 transition (u> = a 2 {\ — \) in dimen- 
sionless units) in the presence of the vacuum state. The smeared two-point function, S(t — s), of equation 
(25) will be employed. After integrating out the angle variables, we get 

^ f°° P d P -i P (t-») 



3tt 2 J [(p/a) 2 + 9/4] 4 

From this point, the investigation of the spontaneous emission process, in particular of the type of decay- 
curve, and of the life-time of the 2P state, can be addressed with the help of mathematical methods for 
Volterra equations [10]. 



4 Conclusions and outlook 

We have presented a simple model of the spontaneous emission of light. While the restriction to the one- 
excitation (photon) transitions is the most obvious limitation of this model, it offers, on the other hand, 
a generalization with respect to the allowed initial states of the quantum electromagnetic field. It asserts 
that the characteristics of the emission process are not intrinsic or unique properties of atoms, but rather 
they depend on the initial state of the quantum electromagnetic field 22 . Two future lines of research seem 
interesting: firstly one might ask, whether the emission process in the distant (cosmological) past was 
similar to the one we observe in laboratories nowadays. The two-point function of the initial state of the 
quantum electromagnetic field, at that epoch, was certainly different from the (Minkowski-) vacuum two- 
point function. The second direction, which we regard as promising, is to endeavor to use simple quantum 

This point is obvious for example from the explicit results of [8]. 
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systems (e.g. some probe-atoms) as detectors of subtle non-classical properties of quantum electromagnetic 
fields. The temporal dependence of the excited-state's amplitude, if measured with sufficient resolution, 
might reveal fingerprints of the, for instance, negative energy-densities of the field. 

Acknowledgements The authors would like to thank prof. K. Fredenhagen for his remarks on the structure of the 
model, as well as for drawing our attention to the works of V. Bach. One of the authors (P. M.) acknowledges the 
financial support of DFG and the hospitality of DESY. 

A Dimensionless units 

In this paper we use dimensionless units. They are defined as follows: 

2 ^ 

E = mc Ed, x = — X£>, 

mc 

h 

t= »t D , p = mcp D , 

mc 1 

^M-J ' A=AD y— 

where the subscript D denotes the dimensionless quantities 23 . The only constant which can appear in the 
equations is the fine-structure constant a. If we take the interaction term as in the equation (3), then the 
quantities corresponding to the free electromagnetic field, in particular: the two-point function, should not 
contain a. As the unit of time is particularly important in the context of this paper, we note that to = 1 
corresponds to h/m e c 2 w 1.25 • 10 _21 s. 

References 

[1] V. Weisskopf and E. Wigner, Berechnung der natiirlichen Linienbreite auf Grand der Diracschen Lichtteorie, Z. 
Phys. 63, 54 (1930). 

[2] M.O. Scully and M.S. Zubairy, Quantum Optics (Cambridge University Press, 1997). 
[3] R. Haag, Local Quantum Physics, 2 ed. (Springer, 1996). 

[4] J. Polzik, E.S. Carri and H.J. Kimble, Spectroscopy with squeezed light, Phys. Rev. Lett. 68, 3020 (1992). 

[5] J. Weinacht, T.C. Ahn and P.H. Bucksbaum, Measurement of amplitude and phase of a sculptured Rydberg wave 

packet, Phys. Rev. Lett. 80, 5508 (1998). 
[6] J.S. Bell and J.M. Leinaas, Electrons as accelerated thermometers, Nucl. Phys. B 212, 131 (1983). 
[7] J. Bach, V. Frohlich and I.M. Sigal, Spectral analysis for systems of atoms and molecules coupled to the quantized 

radiation field, Commun. Math. Phys. 207, 249 (1999). 

As in the paper we use only the dimensionless expressions, this subscript is dropped for brevity. 

Copyright line will be provided by the publisher 



15 



[8] I.V. Bondarev and Ph. Lambin, Spontaneous decay dynamics in atomically doped carbon nanoutubes, Phys. 
Rev. B 70, 035407 (2004). 

[9] M. Raimond, J.M. Brune and S. Haroche, Colloquium: Manipulating quantum entanglement with atoms and 

photons in a cavity, Rev. Mod. Phys. 73, 565 (2001). 
[10] P. Linz, Analytical and Numerical Methods for Volterra Equations (SIAM, Philadelphia, 1985). 
[11] P Marecki and N. Szpak, Spontaneous emission of light from atoms: quantitative analysis of the evolution 

equation, in preparation. 



Copyright line will be provided by the publisher 



